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Summary 


r 


This  technical  report  describes  the  results  of  a  6-month  research  project  conducted  by  Dr.  Ou  Ma  and  one 
of  his  graduate  students  at  the  New  Mexico  State  University  under  the  AFRL/AFOSR  grant  #  FA9550- 
06-1-0284  starting  in  summer  of  2006.  In  the  project,  a  methodology  of  on-orbit  identification  of 
spacecraft  inertia  properties  (i.e.,  the  mass,  the  location  of  the  mass  center,  and  the  inertia  tensor)  using  a 
robotic  arm  was  investigated.  The  extension  of  the  method  to  using  a  solar  array  (instead  of  a  robot)  was 
also  preliminarily  investigated.  The  investigation  was  done  in  theory  and  by  simulation  only  because  of 
the  time  and  funding  limitations.  The  ultimate  goal  of  this  research  is  to  enhance  the  control  and  operation 
of  spacecraft  for  more  challenging  autonomous  on-orbit  service  missions  in  the  future. 

The  investigated  method  makes  use  of  an  onboard  robotic  arm  to  change  the  inertia  distribution  of  the 
spacecraft  system.  As  the  result  of  the  inertia  redistribution,  the  velocity  of  the  spacecraft  system  will 
change  correspondingly.  Since  the  velocity  change  is  measurable  and  the  inertia  change  of  the  robotic  arm 
is  precisely  computable,  the  inertia  parameters  of  the  spacecraft  body  become  the  only  unknown  in  the 
momentum  equations  and  hence,  can  be  identified.  To  treat  the  problem  as  a  linear  identification  problem 
(easier  to  handle  and  solve),  the  problem  has  to  be  formulated  in  two  steps  (or  two  sub-problems).  The 
first  step  is  to  identify  the  mass  and  mass  center  of  the  spacecraft;  and  the  second  step  is  to  identify  the 
inertia  tensor  of  the  spacecraft.  The  largest  advantages  of  this  method  are:  1)  it  does  not  consume  fuel 
because  an  onboard  mechanical  subsystem  is  energized  by  solar  power;  2)  it  requires  to  measure  steady- 
state  velocities  only  (requiring  no  acceleration  and  force  measurements)  because  it  is  based  on  the 
momentum  conservation  law;  3)  it  will  not  be  affected  by  any  internal  energy  dissipation,  which  is  very 
difficult  to  predict  otherwise.  On  the  contrary,  the  existing  thrust  force  based  inertia  identification 
methods  will  not  have  these  desirable  advantages. 

Most  of  the  research  results  reported  in  this  document  has  been  orally  presented  to  and  discussed  with  the 
Dynamics  and  Control  Group  of  AFRL/VSSV  at  Kirtland  Air  Force  Base  on  September  27,  2006. 
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1. 


INTRODUCTION 


1,1  Goals  and  Objectives 

This  research  is  aimed  at  developing  a  new  method  for  on-orbit  identifying  inertia  properties  of  spacecraft. 
The  ultimate  goal  is  to  enhance  the  control  and  operation  of  spacecraft  for  more  challenging  autonomous 
on-orbit  service  missions  in  the  future. 

This  document  summarizes  the  results  of  a  short-term  research  project  conducted  by  the  authors  of  the 
document  at  NMSU  in  the  summer  and  fall  of  2006.  The  objectives  of  the  project  were: 

a)  Develop  a  method  of  using  an  onboard  robotic  arm  (or  arms)  to  identify  the  inertia  parameters  of 
spacecraft.  Conduct  a  thorough  literature  survey  about  the  related  techniques. 

b)  Perform  a  theoretical  study  of  the  method  in  fundamental  aspects  such  as  robustness,  accuracy, 
efficiency,  or  any  other  potential  issues.  Develop  corresponding  solution  techniques  should 
issues/difficulties  be  discovered. 

c)  Implement  a  simulation  model  on  ADAMS  and  data  analysis  program  on  Matlab  to  support  the 
research.  The  simulation  model  should  include  a  rigid  spacecraft,  a  6-DOF  rigid-body  robotic  arm, 
and  sensor  models. 

d)  Explore  the  possibilities  of  following  the  same  idea  but  using  other  means  for  inertia  identification. 
For  example,  use  solar  arrays,  antennas,  other  re-configurable  mechanisms,  or  fuel  level,  as 
opposed  to  use  a  robotic  arm. 

e)  Proposed  a  concept  of  experimental  verification  of  the  proposed  method  in  the  lab  and  on  orbit. 

The  concept  should  be  based  on  the  available  lab  infrastructure  at  either  NMSU  or  AFRL. 

These  objectives  will  help  AFRL  understand  the  proposed  new  method  and  its  potential  benefits  and 
limitations  with  respect  to  the  Air  Force  missions  in  space. 


1.2  Background  and  Application 

Spacecraft  inertia  parameters  can  change  in  orbit  for  many  reasons  such  as  fuel  consumption,  hardware 
reconfiguration,  payload  deployment,  capturing  of  a  flyer,  docking  with  a  satellite,  or  some  mechanical 
malfunctions  like  an  unexpected  deployment  failure.  Spacecraft  state-estimator  or  control  system  needs  to 
know  the  correct  inertia  parameters.  For  example,  after  a  servicing  spacecraft  docked  to  a  target  satellite 
for  servicing,  the  active  vehicle  has  to  stabilize  the  compound  two-vehicle  system  before  a  service  job  can 
be  performed.  Such  a  stabilization  operation  cannot  be  done  optimally  (in  the  sense  of  fuel  economy,  time 
efficiency,  dynamic  impact,  etc.)  if  the  active  vehicle  does  not  have  knowledge  of  the  new  inertia 
properties  of  the  compound  system. 

Robotics-based  (unmanned)  satellite  on-orbit  servicing  has  been  gaining  increasing  attention  in  the 
international  space  community.  Several  technology  development  and  demonstration  missions  have  been 
either  done  recently  or  planned  in  the  near  future  around  the  world  [1-6].  There  are  significant  advantages 
for  having  a  robotics-based  mission  to  rescue  or  service  a  satellite  in  space  for  economical  or  emergency 
purposes.  However,  enabling  technologies  required  for  autonomous  on-orbit  servicing  have  not  been  all 
mature.  Many  research  activities  are  still  going  on  around  the  world  currently.  Japan  completed  an  on- 
orbit  technology  demonstration  mission  ETS-7  in  1999  [1].  In  that  mission,  docking  and  simple  robotic 
operations  were  tested  under  the  condition  that  all  details  of  spacecraft  are  perfectly  known.  DARPA  and 
USAF  is  currently  developing  a  technology  demonstration  mission  through  the  Orbital  Express  Program 
[2].  In  that  mission,  more  advanced  robotics  and  docking  operations  including  capturing  and  refuelling 
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will  be  exercised  on  orbit.  However,  the  servicing  spacecraft  still  perfectly  knows  the  serviced  satellite 
and  the  relative  motion  between  the  two  spacecraft  is  still  gentle.  Germany  and  Canada  are  also  jointly 
developing  a  spacecraft  rescue-and-service  mission  called  TECSAS  [4-5].  In  that  mission,  the  rescuing 
spacecraft  will  be  launched  two  years  after  the  target  satellite,  which  will  also  have  a  slight  tumbling 
motion  during  the  capture  operation.  In  that  case,  identification  of  target  satellite’s  inertia  properties  may 
become  desirable.  Future  on-orbit  rescue  and  service  missions  will  face  more  aggressive  scenarios  such  as 
to  capture  a  tumbling  satellite  with  unknown  inertia.  In  such  a  mission,  quick  identification  of  the  inertia 
properties  of  the  target  satellite  or  the  post-capture  compound  system  will  be  the  key  for  the  capture  and 
post-capture  operations. 

All  of  the  existing  techniques  for  on-orbit  identification  of  inertia  parameters  rely  on  an  active  attitude 
control  system  using  either  thrusters  or  momentum  wheels  to  perform  identification  manoeuvres.  In  this 
research  a  robotics-based  method  is  proposed  and  investigated.  The  method  uses  a  robotic  arm,  which  is 
usually  available  on  servicing  spacecraft,  to  help  identify  the  inertia  properties.  Taking  advantage  of  the 
well  developed  robotics  technology,  the  method  offers  several  benefits  over  the  existing  methods,  as 
described  in  Section  2.8.  The  largest  advantage  of  using  a  robotic  arm  to  do  the  inertia  identification  is 
that  it  requires  no  fuel  consumption  because  an  onboard  robot  is  energized  by  solar  power  as  opposed  to 
fuel.  The  method  can  also  be  extended  to  using  solar  arrays  or  other  motion-controllable  onboard 
mechanisms  (instead  of  a  robot)  to  identify  inertia  properties  of  spacecraft. 


1.3  Literature  Survey 

There  exist  several  methods  for  on-orbit  identification  of  inertia  properties  [7-13].  Although  they  are 
different  in  details,  most  of  them  (i.e.,  [7-12])  use  the  same  fundamental  approach,  which  is  based  on  the 
Newton-Euler  equation  of  motion,  namely, 

'”,v-/»JpJxd»  +  (p5x«))x«)]  =  Sf<  (1) 

Ijd)  +  co  x  Ijto  +  ms ps  x  v  =  It  +  S(p,xf.)  (2) 

where  scalar  ms ,  vector  p  s  and  matrix  I  v  are  the  mass,  position  vector  of  mass  center,  and  inertia  matrix 
of  the  spacecraft  with  respect  to  the  spacecraft  frame,  respectively.  These  three  terms  are  exactly  the 
inertia  properties  to  be  identified.  Other  terms  in  eqs.(l)  and  (2)  are:  v  being  the  acceleration  vector  of  the 
origin  of  the  spacecraft  frame;  co  and  cb  being  the  angular  velocity  and  acceleration  vectors  of  the 
spacecraft  which  have  to  be  measured;  Ep,  x  f,  being  the  resultant  torque  of  the  individual  thrust  forces 
f, ;  and  Zt  being  the  resultant  moment  from  other  external  sources  such  as  gravity-gradient,  residual  air 
drag,  and  solar  pressure.  In  order  to  estimate  the  unknown  inertia  properties  from  equations  (1)  and  (2), 
one  can  convert  the  two  equations  into  a  regression  form  and  then  applies  the  least-square  or  other 

filtering  techniques  to  solve  it  for  the  unknown  components  of  ms  ,  ps  andl  v  (or  I^-1  instead).  It  is  easy 
to  understand  from  the  two  equations  that  such  an  approach  requires  measuring  not  only  the  velocities  but 
also  the  accelerations  and  the  fired  thrust  forces  of  the  spacecraft.  The  procedure  can  be  roughly 
illustrated  in  Fig.l 


Fig.l  Procedure  of  force-based  approach  for  on-orbit  identification  of  inertia  properties 


Precisely  measuring  the  time  history  of  the  magnitude  and  direction  of  each  thrust  force  is  difficult. 
References  [7-10]  basically  followed  the  above  outlined  approach.  References  [1 1]  and  [12]  smartly 
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improved  the  above-mentioned  force-based  approach  by  eliminating  the  nonlinear  term  to  x  Ito  from  the 
regression  equation  by  pre-multiplying  to  x  to  each  term  of  eq.(2).  This  means  that,  instead  of  using  the 
vector  equation,  a  specific  projection  perpendicular  to  to  is  used  to  form  the  regression  system  for 
identification.  This  improves  the  estimation  process  but  the  method  still  requires  known  external 
excitations  (such  as  firing  thrusts).  Without  external  excitations,  the  method  works  only  for  rotating 
spacecraft  with  an  assumption  that  the  initial  rotational  kinetic  energy  is  already  known.  The  work 
reported  in  [12]  also  considered  the  effect  of  the  gravity- gradient  torques  to  the  identification  results. 

The  only  published  work  which  is  fundamentally  different  from  the  above-mentioned  force-based 
approach  is  the  one  reported  by  [13].  That  approach  is  based  on  the  principle  of  conservation  of  angular 
momentum  and  thus,  is  a  momentum-based  approach.  The  method  uses  momentum  wheels  and  the 
associated  attitude  control  system  to  perform  a  set  of  attitude  manoeuvres  and  then  measures  the 
corresponding  attitude,  angular  velocity,  and  wheel  momentum  for  the  identification  process.  Since  the 
motion  (or  degrees  of  freedom)  of  a  momentum  wheel  is  limited  and  the  method  is  based  on  angular 
momentum  equation  only,  it  cannot  identify  the  mass  and  position  of  mass  center  of  the  spacecraft. 

Identification  of  the  inertia  properties  of  spacecraft  using  an  onboard  robotic  arm  is  a  new  approach.  The 
closest  works  published  in  the  literature  would  be  [14-15]  where  the  authors  proposed  methods  of 
identifying  the  inertia  properties  of  a  payload  handled  by  a  robotic  arm  in  space.  The  references  discussed 
both  momentum-based  method  and  force-based  methods  for  identification  of  the  inertia  parameters  of  the 
robot’s  payload. 

The  possibility  of  identifying  spacecraft  inertia  parameters  using  solar  arrays  has  never  been  studied  and 
reported  in  the  open  literature.  Hence,  no  survey  results  in  this  regard  can  be  reported. 
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2.  METHODOLOGY 


2.1  Basic  Concept 

This  subsection  gives  a  brief  idea  about  the  proposed  method  for  those  who  do  not  want  to  go  through  the 
full  mathematical  description  of  the  method.  The  complete  and  detailed  description  is  presented  in  the 
next  few  subsections. 

The  proposed  method  is  based  on  the  principle  of  impulse  and  momentum.  For  easy  understanding  of  the 
concept,  let  us  first  assume  no  external  forces  and  torques  applied  to  the  spacecraft  system  (which 
consists  of  both  the  spacecraft  body  and  the  robotic  arm).  Because  of  no  impulse  from  external  forces  and 
torques,  the  angular  momentum  of  the  system  will  be  conserved,  namely, 

he  =(Is  =  constant  vector  (3) 

where  hc  is  the  angular  momentum  of  the  system;  to  is  the  angular  velocity  of  the  system;  I s  and  I R 
are  the  inertia  matrices  of  the  spacecraft  body  and  the  robotic  arm,  respectively. 


Fig.2  The  arm  moves  to  different  positions  (configurations) 


The  angular  momentum  will  be  conserved  even  if  the  robotic  arm  moves  from  one  configuration  to 
another  with  respect  to  the  spacecraft  body,  as  shown  in  Fig.2.  In  other  word,  at  different  arm 
configurations,  we  have 

hc  =(ls  +Ifli)®i  =(J5  +I*2)(02  =••”  =  (Is  +!*,)«>«  =  constant  vector  (4) 

where  I  ,  I  > ' "  >  I  An  an^  w>  i ,  c#  j , — ,  <a  „  are  the  inertia  matrices  of  the  robotic  arm  and  angular 
velocities  of  the  system  corresponding  to  the  arm  configurations  1,2,  ?  ,  n,  respectively. 

The  inertia  matrix  of  the  robot  at  each  configuration  can  be  computed  using  the  robot’s  dynamics  model 
and  the  known  joint  positions.  The  angular  velocity  of  the  spacecraft  system  can  be  measured  using  rate 
sensors.  Thus,  the  only  unknown  remaining  in  equation  (4)  is  the  inertia  matrix  of  the  spacecraft  body  I  s 
and,  therefore,  it  can  be  solved  from  the  equation. 

The  method  makes  use  of  an  onboard  robotic  arm  to  excite  the  dynamics  of  the  spacecraft  and  then 
measure  the  resulting  angular  velocity  changes  as  input  for  identification  procedure,  which  can  be  briefly 
illustrated  in  Fig. 3. 


Fig.3  Procedure  of  the  new  approach  for  on-orbit  identification  of  inertia  property 
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When  external  forces  and/or  torques  are  not  negligible,  the  resultant  impulse  generated  by  the  resultant 
external  torque  must  be  added  to  equation  (4).  The  angular  impulse  between  two  measurements  can  be 
computed  as 

Ahc  =  |,+l  [St  +  E(p,  x  f,  )}dt  (5) 

where  tt  and  tM  are  respectively  the  times  when  the  z'th  and  (i+l)\h  measurements  were  taken  place. 

2.2  Definitions 

To  facilitate  the  development  and  discussion  of  the  formulation,  we  define  some  basic  terminologies  and 
coordinate  systems  as  follows  (referring  to  Figs.  4  and  5): 

Robot  or  robotic  arm :  the  physical  arm  staring  from  the  shoulder  joint  to  the  end-effector  (i.e.,  the  tip). 
Spacecraft  or  Spacecraft  body,  the  physical  spacecraft  excluding  the  robotic  arm. 

Spacecraft  system :  the  entire  spacecraft  system  including  both  the  spacecraft  body  and  the  robotic  arm. 
Cs  :  mass  center  of  the  spacecraft  body. 

CR :  mass  center  of  the  robotic  arm. 

C  :  mass  center  of  the  spacecraft  system  which  is  also  used  as  the  origin  of  the  orbital  frame. 

Os  :  origin  of  the  spacecraft  frame  which  is  a  fixed  point  on  the  spacecraft  body. 

Fq  :  inertial  frame. 

Fe  :  Earth  frame  which  is  an  inertial  frame  originated  at  the  center  0E  of  the  earth. 

Fs  :  spacecraft  frame  which  is  fixed  to  the  spacecraft  body  and  originated  at  point  Os  of  the 
spacecraft  body. 

ms :  mass  of  the  spacecraft  body. 
mR :  mass  of  the  robotic  arm. 

1 5 :  inertia  matrix  of  the  spacecraft  body  with  respect  to  its  own  mass  center  Cs  .  It  is  a  3D3  constant 
matrix  in  the  spacecraft  body  frame. 

I R :  inertia  matrix  of  the  robotic  arm  with  respect  to  its  own  mass  center  C R .  It  is  a  function  of  the 
configuration  of  the  robot  with  respect  to  the  spacecraft  body. 

ls,c :  inertia  matrix  of  the  spacecraft  body  with  respect  to  the  mass  center  of  the  spacecraft  system  C. 

I  up  '■  inertia  matrix  of  the  robot  with  respect  to  the  mass  center  of  the  spacecraft  system  C. 

Ic :  inertia  matrix  of  the  spacecraft  system  respect  to  its  the  mass  center  C. 

p5 :  position  vector  of  the  mass  center  of  the  spacecraft  body  with  respect  to  the  origin  of  the 
spacecraft  frame.  It  is  the  vector  from  Os  to  Cs . 

p R :  position  vector  of  the  mass  center  of  the  robot  with  respect  to  the  origin  of  the  spacecraft  frame.  It 
is  the  vector  from  Os  to  CR . 

pc :  position  vector  of  the  mass  center  of  the  spacecraft  system  with  respect  to  the  origin  of  the 
spacecraft  frame.  It  is  the  vector  from  Os  to  C  . 
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cs  :  position  vector  of  the  mass  center  of  the  spacecraft  body  with  respect  to  the  mass  center  of  the 
spacecraft  system.  It  is  the  vector  from  C  to  Cs . 

CR :  position  vector  of  the  mass  center  of  the  robot  with  respect  to  the  mass  center  of  the  spacecraft 
system.  It  is  the  vector  from  C  to  Cs  . 

Pf;(0)  is  the  initial  linear  momentum  of  the  spacecraft  system; 
hc(0)  is  the  initial  angular  momentum  of  the  spacecraft  system; 

p c  is  the  linear  momentum  of  the  spacecraft  system  with  respect  to  the  mass  center  C  of  the  system; 

h 

c  is  the  angular  momentum  of  the  spacecraft  system  with  respect  to  the  mass  center  C  of  the  system; 
cs  :  velocity  vector  of  the  mass  center  of  the  spacecraft  body  with  respect  to  the  orbital  frame. 

CR :  velocity  vector  of  the  mass  center  of  the  robot  with  respect  to  the  orbital  frame. 
vs  is  the  linear  velocity  of  the  origin  Os  of  the  spacecraft  frame  Fs  . 

as  is  the  angular  velocity  of  the  spacecraft  body.  This  is  also  the  angular  velocity  of  the  robotic  arm 
when  it  is  locked  in  a  specific  configuration. 

yT 

'  is  the  resultant  external  torque  exerted  on  the  spacecraft  system;  and 

Yf 

i  is  the  resultant  external  forces  exerted  on  the  spacecraft  system. 


Fph 


Fig.4  Illustration  of  the  coordinate  system  for  dynamics  modeling 
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Fig.5  Mass  centers  of  the  robotic  arm,  the  spacecraft  body,  and  the  entire  spacecraft  system 


2.3  Basic  Kinematics  and  Dynamics  Relations 

Based  on  the  definition  of  mass  center  and  referring  to  Fig.5,  we  have  the  following  identities: 

.  mo 

mSCS  +mRCR  =°  0r  CR  = — - CS 

mR 

msps+mRpK 

Pc  _ 

ms  +  mR 

from  which,  we  can  further  derive 


cs  ~  Ps  Pc  -  (Ps  Pr) 


Oh 

S  +  "‘R 


ms  +m 
m 


ms  +mR 


-(Pr-Ps) 


CR  Pr  Pc- 

Ic  =  Is  +  ms  [(cscs)l  ~CSCS  ]+  Ir  +  ml<  [(CrCji)I  —  crcr] 


YYl  772  (  \ 

=  IS  +IR  + — — —  l(Ps  -PR)r(pR  -Pr)1-(Pr -PrXPs  -Pr)7  j 
ms  +mR 

If  the  robot  is  locked  in  a  configuration  without  relative  motion  with  respect  to  the  spacecraft,  then 


vc  =v5  +o>5xpc 
vcs  =c5  =  +to5xp5 

VC«  =Cr  =ViS+(05xp;? 


(6) 


(7) 


(8) 
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The  dynamics  equation  of  the  spacecraft  system 


^-^-ZT  +  Z(P,xfl) 


dt 


dt 


(9) 


or  in  integration  form  as 

hc  =  Icto  =  hc(0)+  |[£t  +  S(p,  xf \)\dt  (10) 

The  following  is  an  important  identity  which  has  been  repeatedly  used  in  the  derivation  of  the  formulation: 

p  x  (to  x  p)  =  (p^pjco  -  (prco)p  =  [(p^p)  -  pp7"  ]co  (11) 

where  1  is  the  3x3  identity  matrix. 


2.4  Method  by  Measuring  both  Linear  and  Angular  Velocities 

This  method  is  derived  based  upon  the  following  assumptions: 

1)  All  the  external  forces  and  torques  exerted  on  the  spacecraft  system  are  ignored. 

2)  Both  linear  and  angular  velocity  of  the  spacecraft  body  is  known  (measured  or  derived  from  other 
measured  data). 

3)  The  configuration  of  the  robotic  arm  is  known  (measured  or  derived  from  measured  data). 

4)  Each  measurement  is  done  when  the  robotic  arm  has  locked  to  a  known  configuration. 


Since  no  external  forces  and  torques,  the  linear  momentum  and  angular  momentum  are  conserved,  namely 


Pc  =Pc(°) 
hc  =  hc(0) 


(12) 


where  Pc(0)  and  h<~(0)  are  the  initial  linear  and  angular  momentums  of  the  spacecraft  system, 
respectively.  Equation  (12)  can  be  expressed  in  details  as 

Pc  =Ps,C+P«,C  =  ™SiS+mR*R  =fns(yS+(aSXPs)  +  mR^S+(aSXPR)  =  Pc(°) 


hC  =hS,C+h^C  =  lS<a  +  cSxmSiS  +1R°>S  +cRxmR*R 

=  I50)5  +c  sxfflj(vs  +  co<j  xps)  +  I*©s  +c  RxmR(ys  +(osxpR) 

_  mHmR  ,  .  ,  . 

=  I5o>5  +- — * - (p5  -P*)x(vs  +to5  xp5) 

ms  +mR 

msmR 

+  - (PS  -P«)X(V5  +«5xPt?) 

ms  +mR 


=  +  -^^-(Ps  - Ptf)x  k  x (Ps  - P*)] 

ms+mR 

mrffln  (  j 

~  \(Ps_Pfl)  (Ps  _  P/?)l  "  (Ps  _  PttXPs 

ms+mR 

=  Ic«s  =  hc(0) 


(13) 


(14) 
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Note  that  the  unknown  inertia  variables  ms  and  p<>  appear  nonlinearly  in  both  equations.  This  is  not 
suitable  for  applying  linear  identification  techniques.  We  can  overcome  the  nonlinear  problem  by  solving 
the  problem  in  two  separate  steps.  The  first  step  is  to  identify  mass  ms  and  position  vector  of  mass  center 
and  the  second  step  is  for  the  inertia  matrix.  Details  are  described  next. 


Step  1:  identification  of  mass  ms  and  position  of  mass  center  p5 
Rewrite  equation  (13)  as 


m. 


me 


vs+<osxpft- 


Pc(°) 


m 


+  »sxps 


=  —V  t 


R  J 


which  is  linear  in  terms  of  the  unknown  terms  mR!  ms  and  p5  . 


(15) 


By  moving  the  arm  into  a  sequence  of  different  configurations,  C^Cj,  ••  -,Cm ,  and  measuring  the 
resulting  linear  and  angular  velocities  of  the  spacecraft  body,  (vs,,G)Sl),  (vSl,to5l),  •  ■  (vSm,toSm) ,  we 

can  construct  a  large  system  of  linear  regression  equations  in  terms  of  the  unknowns  mR  /  ms  and  p5  . 
Then,  solve  these  two  unknowns  from  the  linear  regression  equation  using  least-square  method.  Since 
mR  is  a  known  and  constant,  we  can  find  ms  easily  after  have  found  the  unknown  term  mR  /  ms .  We 
have  two  options  of  writing  the  linear  regression  equation  depending  on  whether  or  not  the  initial  linear 
momentum  pc(0)  is  known.  The  two  options  will  be  described  next: 


Option  1 :  the  initial  linear  momentum  pf  (0)  is  known. 

If  the  initial  momentum  of  the  system  is  known,  then  we  can  construct  the  linear  regression  equations  as 
follows 


where 


Ax  =  b 


(16) 


(17a) 


II 

< 

< 

*Sk 

+  (aSk  X  P  Rk 

Pc  (0)1 

£lk 

LI 

mR  J 

_ 

b*  =- 

0 

0}Sk  (3)  G>sk 

(2)1 

hi 

G 

= 

®s*(3) 

0 

aSk  0) 

-  a  ski2) 

wSk  (0 

0 

k  =  l. 

2,-, 

m 

(17b) 
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where  v5i.  and  CO Sk  are  linear  and  angular  velocities  of  the  spacecraft  body  measured  at  the  Mi 

configuration  of  the  robotic  arm.  Vector  pRk  is  the  position  of  the  mass  center  of  the  entire  robot  when  it 
is  at  the  Mi  configuration,  which  can  be  computed  (see  Section  2.6). 


Option  2:  the  initial  linear  momentum  pc  (0)  is  unknown. 

In  this  case,  we  need  to  eliminate  the  initial  momentum.  This  can  be  done  by  re-writing  equation  (13)  as  a 
momentum  increment  equation  as  follows 

Pc  -Pc(°)  =  w.s(Avst:  +  A<oa  xps)  +  mR(A\Sk  +Acos*  xptt  +to50xApw)  =  0  (18) 


A vSk ,  A(0Sa  and  Ap^  are  the  changes  of  the  linear  velocity,  angular  velocity,  and  mass  center  of  the 
spacecraft  body  from  the  initial  configuration  to  the  Mr  configuration  of  the  robot.  C is  the  initial 

angular  velocity  of  the  spacecraft.  Notice  that  the  resulting  equation  no  longer  has  the  initial  momentum 
term,  from  which,  we  can  form  the  following  new  linear  regression  equation: 


where 


and 


Ax  =  b 


mR/ms 
.  Ps  . 
'A,' 

:  •  > 

A„_ 


A/t  -  [(Av  Sk  +  Ac)  sk  x  pRk 

+  cos0  x  Apw) 

ft*] 

bA  =  -Ava 

0 

-A®a(3) 

A COsk  (2) 

£lk  =  AcojS(fex  = 

A®sit(3) 

0 

-Afl>a(l) 

% 

§ 

<1 

1 

_ i 

A®st(l) 

0 

A U>Sk  -  01 Sk  -e>S0 

Ap  Rk  ~  P  Rk  ~Pro 
k  =  1,  2,  •••,  m. 


(19) 


(20a) 


(20b) 


Note  that  using  the  relative  velocities  AvS/,  and  Ato Sk  for  identification  has  two  advantages  in 

comparison  of  using  the  absolute  velocities  \$k  an<^  *°Sk  '■  0)  relative  velocities  can  be  easily 

obtained  from  integrating  acceleration  data;  (2)  relative  velocities  have  less  measurement  bias  (i.e., 
systematic  errors). 


Step  2:  identification  of  inertia  matrix  Is 

After  completion  of  Step  1,  the  inertia  variables  ms  and  ps  become  known.  Then  we  can  identify  the 

inertia  matrix  Is .  Again,  these  are  two  options  depending  on  whether  the  initial  angular  momentum  is 
known  or  not. 
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Option  1 :  the  initial  angular  momentum  hc(0)  is  known 

Based  on  equation  (14),  we  can  construct  the  following  linear  regression  equation: 

By  =  c 


where 


40,i) 

4(2,2) 

4(3,3) 

4(1,2) 

4(1,3) 

.4(2,3). 


B 


B, 

i 

n 

_ i 

III 

1 - 

2 

i 

S 

CJ 

_ i 

(21) 


(22a) 


B,  = 


a>sk(  1)  0  0  aSk  (2)  ask  (3)  0 

0  <ySjt(2)  0  ®st(l)  0  aski  3) 

0  0  cosk(S)  0  oa(l)  cosk{  2) 


c*  =  hc(0)-IM«>a  — mRms  (Ps-pM)xh5Ax(ps  -pM)] 
ms  +  mR 

k-  1,2,  -,m. 


(22b) 


where  (0Sk  is  the  angular  velocity  of  the  spacecraft  body  measured  at  the  Mi  configuration  of  the  robotic 

arm.  Terms  pRk  and  I Rk  are  the  position  of  the  mass  center  and  the  centroid  inertia  matrix  of  the  entire 
robot  when  it  is  at  the  Mi  configuration,  which  can  be  computed  from  the  robot  model  (see  Section  2.6). 
Term  Is  (/',  j)  is  the  (/j)th  component  of  the  inertia  matrix  Is . 

Option  2:  the  initial  angular  momentum  hc  (0)  is  unknown 

In  this  case,  we  need  to  eliminate  the  initial  momentum.  This  can  be  done  by  re-writing  equation  (14)  as  a 
momentum  increment  equation  as  follows 


he  -hc(0)  -  I5Ao>a 
+ 


4a  + — ((Ps  -P«a)7(Px  -Paa)1~(Ps  "Pm-XPs  “P^a)7) 

ms  +  mR 


A  to 


Sk 


+ 


4a  -4o  + - ~  ~  ((Ps  _P«A)r(Ps  _Pm)1  _(Ps  ~Pt?aXPs  —  P/?A  )^ ) 

ms+mR 

~~ — R  ^  '((Ps  _P«o)r(Ps  _P^o)l_(Ps  “PmXPs  ~ P ro)T)® so 
ms+mR 

0 


co 


50 


(23) 
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where  Aco^  =  co^  —  (%Q  is  the  increment  of  the  angular  velocity  of  the  spacecraft  body  from  the  initial 

configuration  to  the  Ath  configuration  of  the  robotic  arm.  Hence,  we  can  form  the  following  new  linear 
regression  equation: 

By  =  c  (24) 

where 

'4  0,0' 

Is  (2,2) 
s  Is( 3.3) 
y  40,2)  ’ 

4(1,3) 

.4(2,3), 


2.5  Method  by  Measuring  Angular  Velocity  Only 

The  approach  introduced  in  Section  2.4  requires  that  both  the  linear  and  angular  velocities  of  the 
spacecraft  are  known.  Measurement  of  the  linear  velocity  is  not  as  easy  as  the  measurement  of  angular 
velocity.  Usually,  linear  velocity  is  integrated  from  accelerometer  data,  which  will  have  drifting  errors.  In 
this  section,  we  introduce  an  alternative  approach  which  does  not  require  the  known  linear  velocity  of  the 
spacecraft.  This  method  is  derived  based  upon  the  following  assumptions: 

1)  All  the  external  forces  and  torques  exerted  on  the  spacecraft  system  are  ignored. 

2)  The  angular  velocity  of  the  spacecraft  body  is  known  (measured  or  derived  from  measured  data). 

3)  The  configuration  of  the  robotic  arm  is  known  (measured  or  derived  from  measured  data). 

4)  Each  measurement  is  done  when  the  robotic  arm  has  locked  to  a  known  configuration. 

Since  no  external  forces  and  torques,  the  angular  momentum  is  conserved,  namely, 

lc«,s=hc(0)  (26) 

Notice  that  the  inertia  matrix  of  the  spacecraft  system  is  the  sum  of  the  inertia  matrices  of  the  spacecraft 
body  and  the  robotic  arm,  namely, 

Ic  =ls,c  +  lR,C  (27) 

Based  on  the  definitions  of  inertia  matrix  and  mass  center,  we  have 
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4,C  -  1 5  +  mS  [(CSCS  )1  _  CSC S  ] 

(2  \ 

- - - Y  [CP 5  “Pi?)r(Ps  -Ptf)l~(Ps  "PtfXPs  _P^)r] 

(ms+mR)  ) 

4,C  “  4  +  mR  [(cfic«  )1  -  cRcTR  ] 


—  I y?  +mR  - —  2"|[(Ps  _P«)r(PS  “P^)l“(Ps  “P/jXPs  _P^)r] 

U ms+mR )  ) 

where  1  is  the  3?3  identity  matrix.  Substituting  (27)  and  (28)  into  equation  (26),  we  get 
2  2  "I 

4©s  +  mS,mR  +ms”tn  [(p&.  -pfi)r(ps  -Pk)1-(Px  ~PflXPs  -Pfi)rK  =  hc(°)-I««4  ( 

(ms+mR) 

This  is  the  fundamental  equation  for  spacecraft  inertia  identification.  In  the  equation,  mass  ms ,  mass 
center  ,  and  inertia  matrix  I5  are  the  unknowns  to  be  identified  and  all  the  other  variables  can  be 
either  measured  or  computed.  Notice  that  the  inertia  matrix  is  appears  linearly  in  the  equation  but  the 
mass  ms  does  not.  In  order  to  apply  linear  regression  techniques,  all  the  unknown  variables  to  be 
identified  must  be  linear  in  the  regression  equations.  We  can  solve  this  problem  by  defining  a  new 
dimensionless  variable  X  as 

mlmR  +  msm\ 

1=  -  (. 

(ms  +mR) 

Apparently,  X  is  a  function  of  ms .  Since  one  can  always  find  ms  in  terms  of  X  from  (30),  we  will 
identify  X  instead  of  ms .  To  this  end,  equation  (29)  is  re-written  as 

Isws  +AAo>s  =hc(0)-IK<oi.  C 

where 

A  =  [(ps  _ P/j)  (Ps  ~ Pr)1  - (Ps  _ Pi?)(Ps  —  Px)  ]  (--/ 

Since  1 5  is  symmetric,  it  actually  has  only  six  unknown  components.  Thus,  we  can  define  a  minimum  set 
of  unknown  variables  as  follows: 

r  -  ~i 


Ps(  1) 
Ps(  2) 
Ps@) 

4(  U) 
4(2,2) 
4(3,3) 

40.2) 

40.3) 
.4(2,3) 
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(34) 


As  seen  clearly,  equation  (3 1)  is  a  nonlinear  equation  in  terms  of  the  unknown  x,  namely, 

<p(x,  p/;)  =  hc(0)-I^tos 

Obviously  this  is  a  nonlinear  identification  problem  from  which  we  can  solve  for  x  in  terms  of  known 
and  . 


2.6  Inertia  and  Momentum  of  the  Entire  Robotic  Arm 

In  this  section,  we  describe  the  formulation  for  computing  the  inertial  property  (i.e.,  the  mass  center  and 
inertia  matrix)  and  momentum  of  the  entire  robot.  Although  the  mass  center  and  inertia  matrix  of  the 
robot  will  change  with  the  motion  of  the  robot,  they  can  be  precisely  computed  as  long  as  the  inertia 
properties  of  individual  links  of  the  robot  are  known  in  advance.  This  is  because  the  configuration  of  the 
robot  are  fully  controlled  and  known.  First,  the  following  assumptions  are  made: 

1)  The  mass  of  each  body  of  the  robot  is  constant  and  known. 

2)  The  position  of  mass  center  and  inertia  matrix  of  each  body  of  the  robot  are  known.  They  are  constant 
when  they  are  represented  in  coordinate  frame  fixed  to  their  associated  body. 

3)  The  configuration  and  motion  of  the  robot  with  respect  to  the  spacecraft  body  are  known.  This  known 
information  can  be  represented  by  the  joint  angles  and  joint  rates  of  the  robot. 


The  kinematics  of  the  robotic  arm  is  defined  as  shown  in  Fig.6. 


Notation: 

qi  --  Relative  displacements  (state  variable)  of  Joint  i,  expressed  in  Fo,  i  frame  (the  outboard  frame  of  previous  body), 
e,  --  Relative  angular  displacements  of  Joint  I  (represented  by  Euler  angles  or  quaternions),  expressed  in  Fbj.,. 
n  -  Intrabody  vector  of  body  i,  expressed  in  F,  frame  (the  inboard  frame  of  the  i’th  body). 

Ci  --  Position  vector  of  the  center  of  mass  (Q)  of  body  i  with  respect  to  Joint  i,  expressed  in  F,  frame. 


Fig.6  Kinematics  notation  of  the  robotic  arm 
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By  definition,  the  position  of  the  mass  center  and  the  inertia  matrix  of  the  entire  robot  with  respect  to  the 
spacecraft  frame  Fs  can  be  computed  as  follows: 


n 

Xw'P ' 


/= 1 
n 


/=1 


(35) 


(36) 


where  m ,  is  the  mass  of  the  /th  link;  p,  and  I,  are  respectively  the  position  vector  of  the  mass  center  and 
the  inertia  matrix  of  the  /th  body  with  respect  to  the  origin  of  the  spacecraft  frame  Fs  and  also  represented 
in  Fs  .  They  can  be  computed  using  the  following  kinematic  relations 


p,=ZrA+r,,c,  (37) 

k= 1 

I;  =  Ri  *  I( Rf  +  mi  [(Pf  Pi  )1  P/Pf  ] 


for  i  =  1,2,  ■  •  • ,  n.  In  the  above  equations  '  c,  and  'I,  are  the  known  mass  center  and  inertia  matrix  of  the 
/th  body  with  respect  to  the  body  fixed  frame  Ft  and,  therefore,  their  values  will  not  change  during  the 
robot  motion.  Moreover,  vector  '  r,  and  matrix  '  R,  represent  the  relative  position  and  orientation  from 
the  /th  joint  (or  the  /th  body  frame)  to  the  (/+l)th  joint  or  the  tip  of  the  arm,  expressed  in  the  /th  body 
fixed  frame  and  thus,  they  are  constant  too. 

Linear  and  angular  momentums  of  the  entire  robot,  represented  in  the  spacecraft  frame,  can  be  computed 
as  follows: 


n  n 

•r  =ZW'V« 

/=!  /= 1 

/ 

hR  +P'  XW'V«) 

*-«  (38) 


where  x,  and  to,  are  the  absolute  linear  and  angular  velocities  of  the  mass  center  of  the  /th  link  of  the 
robot,  which  can  be  computed  from  the  known  joint  angles  and  rates  using  the  following  recursive 
formulae: 


Vd  =  v„_,  +  wM  x R,_, (M r,_, c,_, )+0,  xR/ c,. 
'0" 
o  \e, 


to,  =  to,.!  +  R, 


for  /  =  1,2 ,•••,« 


(39) 


Note  that  v0  and  to0  are  the  linear  and  angular  velocities  of  the  spacecraft. 


2.7  Technical  Issues  to  be  Addressed 
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The  above-outlined  method  is  still  a  concept.  In  order  to  prove  its  feasibility  for  real  space  application,  it 
has  to  be  thoroughly  studied.  At  this  moment,  the  following  technical  challenges  have  been  identified  for 
the  proposed  research: 

1)  Since  the  observability  of  I  is  closely  dependent  on  the  mechanical  properties  and  the 
configurations  of  the  robotic  arm,  guidelines  need  to  be  established  as  how  to  select  optimal  arm 
poses  such  that  not  only  all  the  inertia  properties  can  be  identified  but  also  the  resulting  numerical 
procedure  be  robust. 

2)  Since  space  robots  are  usually  flexible  (for  minimizing  mass),  the  effect  of  arm  flexibility  onto 
the  identification  method  needs  to  be  looked  into.  This  may  not  be  a  problem  because  the  method 
requires  measuring  steady-state  tas  and  vs  only,  so  that  the  measurement  can  be  done  after  the 
arm  vibration  dies  down. 

3)  The  effectiveness  of  the  method  is  related  to  the  ratio  of  the  arm  inertia  and  the  spacecraft  inertia. 
If  the  arm  is  too  light  to  sufficiently  influence  the  spacecraft,  it  can  hold  a  piece  of  spare  hardware 
(with  known  inertia)  from  the  spacecraft  during  an  identification  maneuver. 

4)  The  issue  of  sensitivity  to  sensor  noise  and  how  to  reduce  it  has  to  be  addressed.  This  is  a 
common  problem  for  any  identification  method.  Among  many  available  techniques,  the  most 
suitable  will  be  identified  and  tailored  to  our  specific  needs. 

5)  The  effect  of  the  gravity-gradient  torque  on  the  proposed  identification  method  needs  to  be 
investigated.  The  gravity-gradient  torque  may  be  non-negligible  if  spacecraft  is  in  a  LEO  orbit.  If 
that  is  the  case,  the  angular  momentum  caused  by  the  gravity-gradient  must  be  considered  in  the 
formulation. 

6)  An  analysis  of  computational  complexity  is  required.  Since  an  inertia  identification  method  is 
usually  a  part  of  a  spacecraft  control  system  (model-based,  adaptive  control),  it  has  to  be  executed 
at  a  fast  rate  with  a  limited  memory  space. 

In  fact,  items  3)  to  5)  are  common  for  any  other  identification  methods  too.  These  problems  and  any 
others,  which  may  be  unveiled  later,  will  be  carefully  investigated.  Based  on  the  investigation,  practical 
solution  measures  will  be  proposed  and  evaluated. 

2.8  Comparison  with  Other  Methods 

From  a  preliminary  study  of  this  new  method,  the  following  pros  and  cons  are  identified: 

-  Since  a  robotic  arm  is  energized  by  solar  energy,  the  method  does  not  consume  any  fuel  of  the 
spacecraft.  On  the  contrary,  the  propulsion  based  methods  consume  fuel. 

-  The  method  is  based  on  the  conservation-of-momentum  principle  and  thus  only  velocity  data  are 
required.  Propulsion  based  methods  require  velocity,  acceleration,  and  force  data  because  they  are 
based  on  equations  of  motion.  Also,  force  and  acceleration  data  are  noisier  than  velocity  data. 

-  The  method  requires  to  measure  steady-state  velocities  only.  This  will  ease  and  make  the 
measurements  much  more  reliable.  Steady-state  measurements  will  have  minimal  errors  caused  by 
the  arm  dynamics  and  spacecraft  structural  dynamics. 

-  Since  the  method  is  based  on  momentum-conservation,  any  unknown  energy  damping  in  the  system 
will  not  affect  the  identification  result.  On  the  contrary,  any  unmodelled  damping  may  corrupt  the 
identification  results  for  a  method  based  on  the  equations  of  motion. 

-  The  method  identifies  the  I  matrix  directly  while  the  force  based  methods  can  only  identify  I”1  or 
I  and  I-1  together,  which  requires  additional  computations. 


19 


The  method  requires  a  robotic  arm  onboard  the  spacecraft.  However,  this  will  likely  be  the  case  for 
spacecraft  designed  for  servicing  a  spacecraft  on  orbit. 

The  method  requires  accurate  inertia  parameters  and  kinematics  of  the  robot.  These  can  be 
guaranteed  before  launch.  After-launch  recalibration  of  robot  parameters  can  be  done  based  on  the 
same  (momentum-conservation  based)  approach  before  the  spacecraft  inertia  properties  have 
changed. 

As  mentioned  in  [15],  the  method  requires  the  spacecraft  to  have  an  initial  angular  momentum  if  the 
inertia  tensor  needs  to  be  identified. 
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3.  SIMULATION  STUDY 


Computer  simulation  is  the  only  means  for  the  research  on  this  short-term  project.  In  this  section  the 
dynamics  simulation  model  and  the  simulation  results  are  described  in  details. 

3.1  Simulation  Model 

The  dynamics  model  for  the  simulation  study  consists  of  a  rigid  spacecraft  and  a  6-DOF  robotic  arm.  The 
kinematics  of  the  system,  represented  by  nine  coordinate  frames,  are  defined  in  Figs.7&8.  The  spacecraft 
frame  Fs  is  fixed  at  the  geometric  center  Os  of  the  spacecraft  body.  The  z'th  local  coordinate  frame  of 
the  robot  is  attached  to  the  /th  link  and  originated  at  the  z'th  joint  of  the  robot  ( z  =  1, 2,  •  •  • ,  6) .  So,  the  frame 
will  move  with  the  z'th  link.  The  tip  frame  F-j  is  attached  at  the  tip  of  the  end-effector,  as  shown  in  Figs. 7 
&8.  The  home  configuration  of  the  robot  (where  all  the  joint  angles  are  zero)  is  shown  in  Fig.7,  namely, 
the  arm  being  completely  folded  to  the  top  surface  of  the  spacecraft  body.  The  attaching  point  of  the  arm 
on  the  spacecraft  body  is  also  shown  in  Fig.8.  The  kinematics  and  dynamics  parameters  of  the 
manipulator  implemented  in  the  simulation  model  are  listed  in  Table  1. 

The  model  has  been  implemented  on  ADAMS  for  dynamic  simulation  and  on  Matlab  for  analysis.  The 
graphics  model  of  the  system  on  ADAMS  is  shown  in  Fig. 9.  Simulation  results  obtained  from  the  model 
are  presented  in  the  next  few  subsections  along  with  discussions  about  the  findings. 


x7  Z5,X6 
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Fig.  8  Kinematic  dimensions  of  the  spacecraft  and  robotic  arm 


Table  1  Kinematics  and  Dynamics  Parameters  of  the  Spacecraft-Arm  System 


Body 

Parameters 

Unit 

Value 

Note 

Spacecraft 

m 

1.0,  1.0,  1.0 

ms 

kg 

1000.0 

SC  5 

m 

[o  0  of 

Sh 

kg-m2 

diag(l  66.67,  166.67,  166.67) 

Link  1 

a\  . 

m 

0.2,0.0,0.025 

m. 

kg 

7.597 

‘e, 

m 

[0  0  a,/2f 
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'I, 

kg-m2 

Computed 

Link  2 

a2  >  ^2  >  $2 

m 

1.0,  0.2,0.025 

m2 

kg 

34.885 

2c2 

m 

2i2 

kg-m2 

Computed 

Link  3 

<*3,b 3,03 

m 

1.0,  0.2,0.025 

m3 

34.885 

3«3 

m 

[flj/2  fc3/ 2  Of 

3I3 

kg-m2 

Computed 

Link  4 

(dummy  link) 

m 

0,  0,0 

mA 

kg 

0 

% 

m 

[0  0  of 

4i4 

kg-m2 

Computed 

Link  5 

(dummy  link) 

m 

0,0,0 

m5 

kg 

0 

5CS 

m 

[0  0  of 

5i5 

kg-m2 

Computed 

Link  6 

(end-effector) 

a6  >  ^6  >  06 

m 

0.5,  0.0,  0.025 

m6 

kg 

15.817 

6c6 

m 

[0  0  a4/2f 

6i6 

kg-m2 

Computed 

Note  that  the  lengths  of  the  4th  and  5th  links  are  zero  because  the  last  three  revolute  joints  are  integrated 
into  a  3-DOF  spherical  joint. 


23 


Fig. 9  Graphics  model  of  the  spacecraft-arm  system  on  ADAMS 


3.2  Effect  of  Arm/Spacecraft  Mass  Ratio 

Since  the  method  is  to  use  a  robotic  arm  to  excite  the  dynamics  of  the  spacecraft  so  that  the  inertia 
properties  can  be  identifiable,  it  is  naturally  concerned  that  the  robotic  arm  has  to  have  sufficiently  large 
inertia  in  order  to  be  able  to  excite  the  dynamics  of  the  spacecraft.  To  address  this  concern,  a  few 
examples  are  illustrated  in  this  subsection  to  see  how  the  identification  results  vary  with  respect  the  ratio 
between  the  arm  mass  and  the  spacecraft  mass. 

In  theory,  the  identification  result  has  nothing  to  do  with  the  arm/spacecraft  mass  ratio.  In  other  words, 
the  unknown  inertia  parameters  can  always  be  100%  accurately  solved  from  the  momentum  equations 
regardless  what  the  mass  ratio  is  as  long  as  the  required  velocity  data  are  perfect.  However,  this  is  not  true 
in  practice  because  the  velocity  data  can  never  be  perfectly  accurate.  In  fact,  the  effect  of  the  mass  ratio  is 
really  in  the  measurement  accuracy  of  the  velocity  data  as  opposed  to  in  the  identification  solution.  In 
general,  velocity  measurement  will  have  larger  relative  error  when  the  mass  ratio  is  small  because  of  the 
correspondingly  less  signal/noise  ratio.  This  larger  error  in  velocity  data  will  in  turn  lead  to  larger 
identification  error.  Based  on  this  argument,  we  will  add  some  errors  in  the  velocity  data  in  order  to 
investigate  the  effect  of  arm/spacecraft  mass  ratio  on  the  performance  of  the  identification  method. 

In  all  the  examples  of  this  subsection,  ten  velocity  measurements  are  taken  at  the  configurations  equally 
separated  along  the  path  defined  in  each  example.  A  bias  error  defined  as  a  percentage  of  the  true  velocity 
value  is  added  to  the  measured  linear  velocity  or  angular  velocity  (but  not  both,  otherwise  the  error  will 
be  cancelled).  The  corresponding  errors  in  the  identified  parameters  are  computed  and  plotted. 

All  the  examples  show  that  the  identification  errors  are  insensitive  to  the  arm/spacecraft  mass  ratio  as 
long  as  the  ratio  is  more  than  0.05.  This  suggests  that  the  mass  of  the  arm  does  not  need  to  be  larger  than 
5%  of  the  spacecraft’s  mass  in  order  to  get  good  identification  quality. 
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Example  3.2,1: 

In  this  case  the  arm  is  commanded  to  go  along  the  general  trajectory  A  defined  in  Section  3.3  (see 
Fig.  13).  Ten  velocity  measurements  are  taken  at  the  configurations  equally  separated  along  the  path.  The 
spacecraft  body  is  initially  spinning  about  the  Z  axis  at  a  rate  of  5  deg/s.  Based  on  the  velocity  data,  the 
inertia  parameters  of  the  spacecraft  are  identified.  The  identification  errors  are  plotted  in  Fig.  10. 


x  10'9 


Fig.  10  Identification  errors  for  different  arm-spacecraft  mass  ratio 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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Example  3.2.2: 

This  example  is  the  same  as  Example  3.2.1  except  Option  2  is  used  to  solve  the  identification  problem  where 
the  initial  momentum  is  assumed  unknown  and  thus,  the  relative  velocities  are  used  in  the  process.  The 
identification  errors  are  plotted  in  Fig.l  1 .  As  we  can  see,  when  the  mass  ratio  is  small,  the  resulting 
identification  errors  by  using  Option  2  is  much  larger  that  using  Option  1 .  This  tells  us  that  we  should  use 
Option  1  to  solve  the  identification  problem  if  condition  allows. 


x  io-4  Method  Option  =  2;  Initial  velocity  =  [8000  0  0  0  0.0011  0.087]  (nVs,  rad/s) 
4 


3 


2 


1 


0 

0  0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5 

Robat/satellite  Mass  Ratio 
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Fig.l  1  Identification  errors  for  different  arm-spacecraft  mass  ratio 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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Example  3.2.3: 


In  this  case  the  arm  is  commanded  to  go  along  the  general  trajectory  B  defined  in  Section  3.3  (see  Fig.  14). 
With  this  arm  maneuver,  an  initial  spinning  about  the  Z  axis  only  will  not  be  enough  to  identify  all  the 
moments  of  inertia.  Thus,  a  tiny  initial  rate  about  the  Y  axis  is  added  to  the  initial  velocity.  The  0.001 1 
deg/s  rate  in  Y  axis  is  equivalent  to  the  angular  rate  of  a  90-min  Earth  orbit.  The  identification  errors  from 
Option  1  are  plotted  in  Fig.  12. 


Fig.  12  Identification  errors  for  different  arm-spacecraft  mass  ratio 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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3.3  Effect  of  Arm  Maneuvers 


The  method  requires  the  robotic  arm  to  move  to  different  configurations  in  order  to  change  the  velocity  of 
the  spacecraft.  It  is  then  understandable  that  these  configurations  should  be  as  different  as  possible  in 
order  to  maximize  the  resulting  new  mass  distribution  and  velocity  variations.  Moreover,  to  be  able  to 
identify  all  the  components  of  the  inertia  tensor,  the  velocity  changes  must  be  in  all  directions.  Therefore, 
the  motion  trajectory  or  path  should  span  as  large  space  as  possible.  In  our  investigation,  four  typical  arm 
motion  trajectories  are  defined  for  simulation  tests.  They  are  described  below. 

Trajectory  A  (seeFig.13) 

It  is  the  most  general  maneuver  having  the  largest  motion-span  or  arm  maneuvering.  The  arm  first 
stretches  all  the  way  upward;  then  it  swings  down  to  the  horizontal  plane  while  keeping  fully  stretched 
out;  then  it  swings  360  degrees  on  horizontally.  The  starting  configuration  of  the  trajectory  is 
0  =  [o  0  0  0  0  O]  and  the  ending  configuration  is  0  =  [l 80°  -180°  180°  -180°  180°  0°]. 

A  few  typical  configurations  of  the  trajectory  are  shown  in  Fig. 13.  The  figure  also  plots  the  time  histories 
of  the  3-D  linear  and  angular  velocities  along  the  trajectory. 

Trajectory  B  (see  Fig.  14) 

This  trajectory  has  the  second  largest  motion-span  or  arm  maneuvering.  The  arm  is  fully  stretched  out  and 
swings  in  parallel  to  the  top  face  of  the  satellite.  The  starting  configuration  of  the  trajectory  is 

0  =  [o  0  180°  -180°  0  o]  and  ending  configuration  is  0  =  [l  80°  0  180°  -180°  0  0  .  The 

starting  and  ending  configurations  are  shown  in  Fig.  14.  The  figure  also  plots  the  time  histories  of  the  3-D 
linear  and  angular  velocities  along  the  trajectory. 

Trajectory  C  (see  Fig.  15)  , 

This  trajectory  has  the  third  largest  motion-span  or  arm  maneuvering.  The  arm  stretches  all  the  way  up 
and  then  folds  all  the  way  down.  The  starting  configuration  of  the  trajectory  is  0  =  [0  0  0  0  0  0] 

and  ending  configuration  is  0  =  [o  -90°  180°  - 180°  0  o].  The  starting  and  ending  configurations 

of  the  arm  are  shown  in  Fig.  15.  The  figure  also  plots  the  time  histories  of  the  3-D  linear  and  angular 
velocities  along  the  trajectory. 

Trajectory  D  (see  Fig.  1 6) 

This  trajectory  has  the  smallest  motion-span  or  arm  maneuvering.  The  arm  is  fully  stretched  upright  and 
spin  by  the  first  joint  (along  the  Z  axis).  The  starting  configuration  of  the  trajectory  is 

0=0  -90°  180°  -180°  0  o]  and  the  ending  configuration  is 

0=180°  -90°  180°  -180°  0  o],  as  shown  in  Fig.  16.  The  figure  also  plots  the  time  histories  of 

the  3-D  linear  and  angular  velocities  along  the  trajectory. 

The  velocity  profiles  associated  with  the  above-mentioned  trajectories  give  a  good  idea  about  how  good 
they  are  for  the  identification.  In  general,  a  trajectory  is  good  if  it  causes  large  velocity  changes  in  all  the 
6  components  when  the  arm  moves  along  the  trajectory,  such  as  trajectory  A.  a  trajectory  is  not  good  if 
one  or  more  velocity  components  are  not  changed  much  along  the  trajectory,  such  as  trajectory  D. 

Among  the  four  typical  trajectories,  trajectory  A  is  the  best  one  for  inertia  identification  because  it  has  the 
largest  motion  span  and  trajectory  D  is  the  worst  one  because  it  has  the  smallest  motion  span.  The  other 
two  trajectories  are  in  the  middle.  This  point  can  be  seen  from  the  example  simulations  presented  in  the 
next  few  subsections. 
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Trajectory  A:  The  arm  moves  very  generally 

starting  at  Initial  configuration:  0  =  [0  0  0  0  0  0]  and 

ending  at  the  configuration  0  =  [l  80°  -180°  180°  -180“  180°  0°] 


"  ^ -j  ~  ^  ^  ®  StoPP6*5 
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Fig.  13  The  initial,  two  middle,  and  the  ending  configurations  of  trajectory  A 
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Trajectory  C:  the  arm  is  fully  stretched  out  and  swings  horizontally 

starting  at  Initial  configuration:  0  =  [o  0  180°  -180°  0  o]  and 

ending  at  the  configuration  0  =  [l  80"  0  180°  -180"  0  o] 
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Fig.  14  The  initial  and  final  configurations  of  trajectory  B 
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Trajectory  C:  the  arm  moves  up  and  down  vertically 

starting  at  Initial  configuration:  0  =  [0  0  0  0  0  0]  and 

ending  at  the  configuration  0  =  [o  -90°  180°  —180°  0  o] 
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Fig.  15  The  initial  and  fully  stretched-up  configurations  of  trajectory  C 
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Fig.  16  The  initial  and  final  configurations  of  trajectory  D 


Trajectory  D:  the  arm  is  fully  stretched  out  and  swings  horizontally 

starting  at  Initial  configuration:  0  =  [o  -90°  180°  -180°  0  o]  and 

ending  at  the  configuration  0  =  [l 80°  -90°  180°  -180°  0  o] 
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3.4  Sensitivity  to  Bias  in  Velocity  Measurements 


Velocity  measurements  always  have  errors  which  will  cause  corresponding  errors  in  the  identification 
results.  To  access  the  sensitivity  to  the  measurement  bias,  a  few  examples  are  illustrated  in  this  section.  In 
all  examples,  ten  velocity  measurements  are  taken  at  the  configurations  equally  separated  along  the  path 
defined  in  each  example.  A  bias  error  (defined  as  percentage  of  the  true  velocity  value)  is  added  to  the 
measured  linear  velocity  or  angular  velocity  (but  not  both).  The  corresponding  errors  in  the  identified 
parameters  are  computed  and  plotted.  The  plotted  identification  errors  are  defined  as  the  differences 
between  the  identified  values  of  the  inertia  parameters  and  the  true  values  of  the  same  inertia  parameters. 

The  plots  show  that  the  errors  in  the  identified  inertia  parameters  caused  by  the  measurement  bias  error 
are  negligible.  The  error  in  the  identified  position  of  mass  center  is  also  very  small.  The  errors  in  the 
moments  and  products  of  inertia  are  less  than  1/5  of  the  measurement  bias  error,  which  are  also 
considered  small.  This  clearly  indicates  that  the  measurement  bias  is  not  a  concern  for  the  proposed 
identification  method. 
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Example  3.4.1: 

The  arm  is  commanded  to  go  along  the  general  trajectory  A  defined  in  Section  3.3  (see  Fig.  13).  The 
spacecraft  body  is  initially  spinning  about  the  Z  axis  at  a  rate  of  1.0  deg/s.  The  identification  errors  are 
plotted  in  Figs.  17- 18. 


x  io’7  Method  Option  =  1;  Initial  velocity  =  [8000  0  0  0  0  0.017]  (m/s,  rad/s) 
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Fig.  1 7  Identification  errors  when  different  bias  errors  are  added  to  the  linear  velocity 

(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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Bias  Level  in  measured  angular  velocity  (%) 


Bias  Level  in  measured  angular  velocity  (%) 
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Fig.  18  Identification  errors  when  different  bias  errors  are  added  to  the  angular  velocity 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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Example  3.4.2: 

In  this  example,  the  ami  is  fully  stretched  out  and  swings  in  the  horizontal  plane  along  the  trajectory  C 
defined  in  Section  3.3  (see  Fig.  15).  The  spacecraft  body  is  initially  rotating  at  a  rate  of  1 .0  deg/s  about  the 
Y  axis  (pitch)  and  5.0  deg/s  about  the  Z  axis  (yaw).  The  identification  errors  are  plotted  in  Fig.  19. 


'7  Method  Option  =  1;  Initial  velocity  =  [8000  0  0  0  0  0.017]  (m/s,  rad/s) 
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Fig.  19  Identification  errors  when  different  bias  errors  are  added  to  the  angular  velocity 
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3.5  Sensitivity  to  Sensor  Noise 

Velocity  measurements  always  have  errors  due  to  sensor  noise.  These  random  errors  will  cause 
corresponding  errors  in  the  identification  results.  To  access  the  sensitivity  to  the  sensor  noise,  a  few 
examples  are  illustrated  in  this  section.  In  all  examples,  ten  velocity  measurements  are  taken  at  the 
configurations  equally  separated  along  the  path  defined  in  each  example.  A  random  error  (with  zero  mean 
and  the  standard  deviation  defined  as  a  percentage  of  the  true  velocity  value)  is  added  to  each  of  the 
measured  velocity  component  and  then  the  corresponding  errors  in  the  identified  parameters  are 
computed. 

Simulation  tests  with  the  four  trajectories  all  show  that  the  random  errors  in  velocity  data  can  cause 
significant  errors  in  identified  inertia  parameters.  For  the  four  different  arm  maneuvering  trajectories,  the 
identification  errors  continuously  increases  from  trajectory  A  to  trajectory  D,  with  A  having  the  minimum 
errors  and  D  having  the  maximum  errors.  This  is  mainly  because,  as  discussed  in  Section  3.2,  trajectory  A 
has  the  largest  arm  maneuvering  and  trajectory  D  has  the  minimum  arm  maneuvering.  If  we  can  control 
the  velocity-measurement  errors  within  3%  level,  the  resulting  identification  errors  will  also  be  more  or 
less  about  3%  level  for  trajectory  A.  However,  the  same  3%  measurement  errors  could  cause  over  200% 
identification  errors  if  using  trajectory  D.  This  is  totally  unacceptable.  Therefore,  we  know  that  selection 
of  arm  maneuvering  trajectory  for  the  identification  process  is  very  important. 

For  the  same  trajectory,  the  errors  in  identified  mass  and  mass  center  seems  to  be  less  than  that  in  the 
inertia  tensor.  In  other  words,  identification  of  inertia  tensor  is  more  sensitive  to  sensor  noise  than  that  of 
mass  and  mass  center. 
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Example  3.5.1: 

The  arm  is  commanded  to  go  along  the  general  trajectory  A  defined  in  Section  3.3  (see  Fig.  13).  The 
spacecraft  body  is  initially  spinning  about  the  Z  axis  at  a  rate  of  1.0  deg/s.  The  identification  errors  are 
plotted  in  Fig.20. 


Method  Option  =  1;  Initial  velocity  =  [00000  0.017]  (rrVs,  radfe) 


Fig.20  Identification  errors  when  different  noises  are  added  to  the  angular  velocity  measurement 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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Example  3.5.2: 


In  this  example,  the  arm  is  fully  stretched  out  and  swings  in  the  horizontal  plane  along  the  trajectory  B 
defined  in  Section  3.3  (see  Fig.  14).  The  spacecraft  body  is  initially  rotating  at  a  rate  of  1.0  deg/s  about 
the  Y  axis  (pitch)  and  5.0  deg/s  about  the  Z  axis  (yaw).  The  identification  errors  are  plotted  in  Fig.2 1 . 
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Fig.2 1  Identification  errors  when  different  noises  are  added  to  the  angular  velocity  measurement 
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ReJ.  Error  in  Inertia  Moments  (%)  Error  in  Mass  Center  (m)  Relative  Error  in  Mass  (%) 


In  this  example,  the  arm  is  fully  stretched  out  and  swings  in  the  horizontal  plane  along  the  trajectory  C 
defined  in  Section  3.3  (see  Fig.  15).  The  spacecraft  body  is  initially  rotating  at  a  rate  of  1.0  deg/s  about  the 
Y  axis  (pitch)  and  5.0  deg/s  about  the  Z  axis  (yaw).  The  identification  errors  are  plotted  in  Fig.22. 


Method  Option  =  1;  Initial  velocity  =  [0000  0.017  0.087]  (m/s,  rad/s) 


Velocity-Measurement  Noise  Level  (%) 
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Fig.22  Identification  errors  when  different  noises  are  added  to  the  angular  velocity  measurement 
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Example  3.5.4: 

In  this  example,  the  arm  is  spinning  when  it  is  fully  stretched  out  and  upright  along  the  trajectory  D 
defined  in  Section  3.3  (see  Fig.  16).  The  spacecraft  body  is  initially  rotating  at  a  rate  of  1.0  deg/s  about  the 
Y  axis  (pitch)  and  5.0  deg/s  about  the  Z  axis  (yaw).  The  identification  errors  are  plotted  in  Fig.23. 


Method  Option  =  1;  Initial  velocity  =  [0  00  0  0.017  0.087]  (m/s,  rad/s) 
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Fig.23  Identification  errors  when  different  noises  are  added  to  the  angular  velocity  measurement 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 
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3.6  Effect  of  Gravity-Gradient  Torque 

Not  investigated  (due  to  time  and  funding  limitations) 


3.7  A  Study  of  the  Case  using  Solar-Array  Mechanisms 

Dr.  Robertson  and  Dr.  Pham  of  AFRL  asked  if  the  robotics-based  identification  method  can  be  extended 
to  more  general  cases  where  other  onboard  mechanisms  (as  opposed  to  just  robotic  arms)  can  be  used  to 
identify  the  inertia  properties.  These  mechanisms  may  include  solar  arrays,  antennas,  momentum  wheels, 
etc.  If  this  question  gets  a  positive  answer,  then  the  impact  of  the  identification  method  will  be  much 
larger  because  the  method  will  no  longer  rely  on  a  robotic  arm.  To  fully  answer  this  question,  an  intensive 
research  is  needed  which  is  beyond  the  scope  of  this  short-term  research  project.  In  this  section,  we 
perform  a  preliminary  study  of  the  possibility  of  using  a  solar  array  for  inertia  identification  just  for 
gaining  some  idea  about  the  outcome  of  this  extended  research. 

For  the  simulation  study,  we  modify  the  system  dynamics  model  defined  in  3.1  by  replacing  the  robotic 
arm  with  two  solar  panels.  The  two  solar  arrays  are  identical  and  symmetrically  attached  to  the  spacecraft 
body,  as  shown  in  Fig.24.  Each  array  is  a  2-DOF  mechanism  for  tracking  the  Sun.  The  mechanical 
parameters  of  the  arrays  are  given  in  Table  2. 


Fig.24  Reference  frames  of  the  spacecraft  with  sun-tracking  solar  arrays  (in  home  configuration) 
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Table  2  Kinematics  and  Dynamics  Parameters  of  the  Spacecraft  with  Solar  Arrays 


Body 

Parameters 

Unit 

Value 

Note 

Spacecraft 

as>  bS’CS 

m 

1.0,  1.0,  1.0 

ms 

kg 

1000.0 

s*s 

m 

[o  0  of 

Sl  s 

kg-m2 

Ixx=Iyy=Izz  =(1/6  )msdj 

Solar 

Array  1 

a,  b,  c 

m 

2.0,  1.0,  0.2 

mA 

kg 

50 

(5%  spacecraft  mass) 

’C1 

m 

[0  flj/2  +  c,  Of 

‘I. 

kg-m2 

~mA  diag(a2  +  b2 ,a2 ,b2) 

Calculated  from  geometry 

Solar 

Array  2 

a,  b,  c 

m 

2.0,  1.0,  0.2 

mA 

kg 

50 

(5%  spacecraft  mass) 

2c2 

m 

[0  -a/2-c  Of 

2h 

kg-m2 

mA  diag {a2  +b2  ,a2  ,b2) 

Calculated  from  geometry 

Simulation  tests  showed  that,  by  using  the  solar  arrays,  all  the  inertia  parameters  of  the  satellite  can  be 
identified  precisely  without  any  problems  if  the  velocities  of  the  satellite  can  be  precisely  measured. 
However,  similarly  as  in  the  cases  of  using  a  robotic  arm,  if  the  velocity  data  has  errors  or  sensor  noise, 
the  identified  inertia  parameters  will  also  have  errors.  Two  simulation  examples  are  included  in  the  next 
few  pages  to  demonstrate  this. 

The  simulations  examples  suggest  that  the  identification  of  the  mass  and  mass-center  location  is  not 
sensitive  to  sensor  noise  but  the  identification  of  inertia  tensor  is  very  sensitive  to  sensor  noise.  Therefore, 
in  order  to  make  the  method  practically  feasible,  we  need  to  either  reduce  the  sensitivity  of  the  method  to 
sensor  noise  or  directly  reduce  the  velocity  measurement  errors. 
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Example  3.7.1: 


In  this  example,  the  two  solar  arrays  are  spinning  about  their  own  axes  by  180  degrees.  In  other  words, 
the  joint  angles  are  from  0  =  [o  o]  to  0  =  [o  1 80°  ]  for  each  solar  panel.  Ten  velocity  measurements  are 
taken  at  the  configurations  equally  separated  along  the  path.  The  spacecraft  body  is  initially  rotating  at  1 .0 
deg/s  about  the  Y  axis  and  also  1.0  deg/s  about  the  Z  axis.  Based  on  the  velocity  data,  the  inertia 
parameters  of  the  spacecraft  are  identified.  The  identification  errors  are  plotted  in  Fig. 25.  As  one  can  see, 
the  errors  in  mass  and  mass  centers  are  negligible  but  these  in  inertia  tensor  are  very  large.  This  is  because 
the  angular  maneuvering  of  the  satellite  caused  by  the  solar  arrays  is  not  large  enough.  Therefore,  we  can 
say  that  the  identification  of  inertia  tensor  is  very  sensitive  to  sensor  noise. 


Velocity-Measurement  Noise  Level  (%) 


Velocity-Measurement  Noise  Level  (%) 


o'- 


Fig.25  The  identification  errors  versus  the  velocity  measurement  noise  level 
(Note:  the  error  in  the  position  of  the  mass  center  is  represented  by  the  absolute  error  because  the  mass 
center  is  located  at  the  origin  and  thus  the  relative  error  is  undefined) 


45 


In  this  example,  the  two  solar  arrays  are  spinning  about  their  own  axes  by  180  degrees.  In  other  words, 
the  joint  angles  are  from  0  =  [o  o]  to  0  =  [30°  1 80°  ]  for  each  solar  panel.  Ten  velocity  measurements 

are  taken  at  the  configurations  equally  separated  along  the  path.  The  spacecraft  body  is  initially  rotating  at 
1.0  deg/s  about  the  Y  axis  and  also  1.0  deg/s  about  the  Z  axis.  Based  on  the  velocity  data,  the  inertia 
parameters  of  the  spacecraft  are  identified.  The  identification  errors  are  plotted  in  Fig.26.  As  one  can  see, 
the  errors  in  mass  and  mass  centers  are  negligible  but  these  in  inertia  tensor  are  very  large.  This  is  because 
the  angular  maneuvering  of  the  satellite  caused  by  the  solar  arrays  is  not  large  enough.  Therefore,  we  can 
say  that  the  identification  of  inertia  tensor  is  very  sensitive  to  sensor  noise. 


Velocity-Measurement  Noise  Level  (%) 


Fig.26  The  identification  errors  versus  the  velocity  measurement  noise  level 


4.  IDEA  FOR  EXPERIMENTAL  VERIFICATION 


Although  an  experimental  study  of  the  proposed  technology  is  beyond  the  cope  of  this  short-term  project, 
an  idea  of  testing  the  technology  in  a  lab  environment  is  presented  in  this  section  for  discussion.  The  idea 
is  illustrated  in  Fig. 27. 


Figure  27  Concept  of  the  a  2-D  experiment  testbed  for  testing  the  proposed  technology 


The  testbed  consists  of  a  granite  table,  an  air-bearing  supported  mobile  frame,  a  spacecraft  mockup,  a 
robotic  arm,  and  a  2-D  motion  sensing  system.  Note  that  the  motion  of  each  robot  link  including  the  end- 
effector  does  not  need  to  be  measured  because  it  can  be  calculated  from  the  well  controlled  joint  motions 
instead.  It  leaves  that  only  the  motion  of  the  spacecraft  mockup  or  the  mobile  platform  needs  to  be 
measured.  Since  this  is  can  be  considered  as  a  rigid-body  motion  in  the  horizontal  plane,  it  is  not  difficult 
to  measure.  We  have  several  options:  1)  use  a  visual  sensor  to  track  the  position  of  the  mobile  frame  and 
then  estimate  the  velocities  by  differentiating  the  position  data  (results  may  be  noisy  and  thus  requiring 
filtering);  2)  use  two  to  three  accelerometers  to  measure  the  accelerations  of  some  specific  points  of  the 
mobile  frame  and  then,  estimate  the  velocity  by  integrating  the  acceleration  data;  3)  use  a  rate  gyro  to 
measure  angular  velocity  directly  but  the  linear  velocity  still  has  to  be  measured  by  a  visual  sensor  or 
accelerometers.  Since  the  negligible  friction  between  the  air  bearings  and  the  granite  table,  the  spacecraft 
and  the  robotic  arm  will  behave  like  free  floating  on  the  top  of  the  granite  table.  This  will  approximate  a 
2-D  (3  degrees  of  freedom)  free  floating  condition  in  the  space. 

In  the  experiment,  we  give  a  slight  push  to  the  mobile  frame  in  order  to  generate  an  initial  motion  (sliding 
and/or  spinning).  Then,  we  command  the  robot  to  go  through  a  series  of  new  configurations  and  measure 
the  linear  and  angular  velocities  at  each  of  these  new  configurations.  If  we  use  visual  position  sensors,  the 
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velocity  can  be  estimated  by  differentiating  the  measured  position  data  (the  results  may  be  noisy).  If  we 
use  accelerometers,  the  velocity  can  be  estimated  by  integrating  the  acceleration  data.  Rate  gyro  can 
direct  measure  angular  velocity. 

Since  this  is  a  2D  experimental  system,  we  can  identify  only  4  out  of  the  1 0  inertia  parameters  of  the 
spacecraft  mockup.  The  identifiable  parameters  are  the  mass,  the  x  andy  coordinates  of  the  mass  center, 
and  the  second  moment  of  inertia  about  the  z  axis.  However,  this  should  be  enough  to  prove  the  principle 
of  the  proposed  identification  method  and  assess  the  feasibility  and  performance  of  the  proposed 
technology. 

An  important  requirement  for  the  testbed  is  that  the  inertia  properties  of  the  robotic  arms  and  the  mobile 
platform  must  be  known  in  advance.  This  is  also  a  condition  of  the  proposed  inertia  identification  method 
in  a  real  application  mission.  This  requirement  will  not  be  difficult  to  meet  because  the  hardware 
components  are  specially  designed,  manufactured,  and  integrated  in  the  lab.  We  can  always  carefully 
measure  and  verity  these  properties  of  the  hardware  in  the  lab. 
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5.  CONCLUSIONS  AND  RECOMMENDATION 

The  following  have  been  accomplished  in  this  short  project: 

a)  Developed  the  method  of  using  an  onboard  robotic  arm  (or  arms)  to  identify  the  inertia  parameters 
of  a  flying  spacecraft.  Two  variations  of  the  method  are  proposed;  one  requires  measuring  both 
linear  and  angular  velocities  and  the  other  requires  measuring  the  angular  velocity  only.  The  latter 
is  a  nonlinear  parameter  identification  problem  and  thus,  much  more  difficult  to  solve.  The  focus  of 
this  project  is  on  the  former  only. 

b)  Implemented  a  simulation  model  on  ADAMS  and  data  analysis  program  on  Matlab  to  support  the 
research.  The  simulation  model  consists  of  a  rigid  spacecraft  and  a  6-DOF  rigid  robotic  arm. 

c)  Performed  a  series  of  simulations  to  test  the  new  method  for  some  basic  concerns  such  as 
robustness,  accuracy,  and  limitations. 

d)  Explored  the  idea  of  using  a  pair  of  solar  arrays  (instead  of  a  robotic  arm)  to  identify  the  same 
inertia  parameters. 

e)  Proposed  a  concept  of  experimental  verification  of  the  proposed  method  in  lab  environment. 

From  the  simulation  tests,  the  following  findings  are  concluded: 

1)  The  method  of  using  an  onboard  robotic  arm  to  identify  the  inertia  properties  of  a  flying  spacecraft 
is  feasible  in  theory  (note:  “in  theory”  because  it  has  not  been  experimentally  verified). 

2)  The  idea  of  using  solar  arrays  (instead  of  a  robotic  arm)  to  identify  the  inertia  properties  of  a 
spacecraft  is  also  possible  in  theory. 

3)  In  order  to  treat  the  problem  as  a  linear  identification  problem  (much  easier  to  solve),  it  has  to  be 
formulated  in  two  steps  (or  two  sub-problems).  The  first  step  is  to  identify  the  mass  and  mass 
center  of  the  spacecraft;  and  the  second  step  is  to  identify  the  inertia  tensor  of  the  spacecraft. 

4)  The  spacecraft  does  not  require  an  initial  spin  (or  rotation)  for  identifying  the  mass  and  mass 
center.  However,  it  does  require  an  initial  spin  (or  rotation)  to  identify  the  inertia  tensor.  This  initial 
angular  velocity  can  be  unknown  (i.e.,  Option  2  of  the  method). 

5)  When  using  a  robotic  arm  for  inertia  identification,  note  the  following: 

-  For  good  identification  performance,  the  arm  mass  should  be  no  less  than  5 — 10%  of  the 
spacecraft  mass. 

-  The  arm  does  not  necessarily  need  6  degrees  of  freedom.  However,  more  degrees  of  freedom 
will  make  the  required  maneuvering  easier. 

-  If  the  spacecraft  initially  rotates  about  one  axis  only,  then  the  arm  motion  should  span  at  least 
two  axes.  If  the  spacecraft  initially  rotates  in  two  axes,  the  arm  motion  can  span  only  one  axis. 

-  The  identification  solution  is  not  sensitive  to  (velocity)  measurement  bias,  e.g.,  10%  bias 
causing  less  than  1%  identification  error. 

-  The  identification  solution  is  sensitive  to  (velocity)  measurement  noise,  e.g.,  1%  noise  can 
cause  over  5%  identification  error. 

-  The  identification  process  is  much  less  sensitive  to  measurement  errors  when  the  initial  angular 
momentum  is  known  (i.e.,  option  1  of  the  method)  than  unknown  (i.e.,  option  2  of  the  method). 

6)  When  using  solar  arrays  for  inertia  identification, 

-  The  solar  arrays  must  be  motion  controllable. 

-  If  the  satellite  initially  rotates  about  only  one  axis,  then  the  solar  arrays  have  to  be  able  to  rotate 
in  two  axes  in  order  to  identify  all  the  inertia  matrix. 
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-  The  method  is  insensitive  to  sensor  noise  for  identifying  mass  and  mass-center  location.  But  it  is 
sensitive  to  sensor  noise  for  identifying  the  inertia  tensor. 

-  Other  findings  of  the  method  should  be  similar  to  these  by  using  a  robotic  arm. 

Recommendation  for  further  study: 

1)  Complete  the  study  of  using  a  solar-array  mechanism 

2)  Continue  the  research  to  address  the  issue  of  sensitivity  to  sensor  noise 

3)  Investigate  the  concern  about  external  disturbances 

4)  Design  an  air-bearing  based  testbed  for  experimental  verification 

5)  Explore  ideas  of  using  other  means,  e.g.,  known  fuel  level 
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6. 
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